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Electrical neurostimulation for chronic pain: 
on selective relay of sensory neural activities in myelinated nerve fibers 
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Abstract —Chronic pain affects about 100 million adults 
in the US. Despite their great need, neuropharmacology and 
neurostimulation therapies for chronic pain have been associ¬ 
ated with suboptimal efficacy and limited long-term success, 
as their mechanisms of action are unclear. Yet current com¬ 
putational models of pain transmission suffer from several 
limitations. In particular, dorsal column models do not include 
the fundamental underlying sensory activity traveling in these 
nerve fibers. We developed a (simple) simulation test bed 
of electrical neurostimulation of myelinated nerve fibers with 
underlying sensory activity. This paper reports our findings 
so far. Interactions between stimulation-evoked and underlying 
activities are mainly due to collisions of action potentials and 
losses of excitability due to the refractory period following an 
action potential. In addition, intuitively, the reliability of sensory 
activity decreases as the stimulation frequency increases. This 
first step opens the door to a better understanding of pain 
transmission and its modulation by neurostimulation therapies. 
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Fig. 1. Pain system builds on a tightly regulated dynamical crosstalk 
between the peripheral nervous system and the brain via the spinal cord. 
Spinal cord stimulation modulates neural activity to suppress chronic pain. 


I. Introduction 

Pain is a protective and adaptive physiological system es¬ 
sential for the survival of many species. However, this system 
is fragile as damage and malfunction of the nervous system 
may divert its function, creating a debilitating disease known 
as chronic pain. Chronic pain affects about 100 million 
adults in the US, with $560-635 billion in annual medical 
expenses and lost productivity [1]. Chronic pain is primarily 
treated with neuropharmacology, which may be inadequate 
or toxic, have negative side effects (for example, addiction 
to narcotics), and lose efficacy after long-term use [2], [3]. 
Alternatively, it is also treated with Spinal Cord Stimulation 
(SCS), an electrical neurostimulation that has potential to 
reduce the need for drugs and produces less negative side 
effects. However, SCS has been associated with suboptimal 
efficacy and limited long-term success as their mechanisms 
of action are unclear [4]. 

Critical to advancing chronic pain treatment is a deeper 
mechanistic understanding of pain transmission and modu¬ 
lation under both normal and pathological conditions, which 
remains largely elusive because the pain system is complex. 
The pain system builds on a tightly regulated dynamical 
crosstalk between the peripheral nervous system and the 
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brain via the spinal cord (see Fig. [T}. Over the past decades, 
detailed computational models have been used to understand 
the effects of electrical neurostimulation on dorsal column 
fibers, a localized but important part of this complex system. 
In particular, efforts have been made (i) to determine the 
response of a nerve fiber to an electrical potential field 
and (ii) to compute the shape of the electrical potential 
field created by stimulation (see [5] for a recent review on 
both research directions). Although these models reproduce 
some observed behaviors, none of these models include the 
fundamental underlying sensory activity (either normal or 
pathological) traveling in these nerve fibers. 

As a first step to address this issue, we developed a 
simulation test bed of extracellular electrical stimulation of 
myelinated nerve fibers in the dorsal column with underlying 
sensory activity. Unlike previous approaches, our approach 
considers the potential interactions in dorsal column fibers 
between neurostimulation-evoked activity and underlying 
sensory activity coming from peripheral nerves. By adding 
this sensory input, we will gain a better understanding of 
different mechanisms for SCS analgesia. 

The paper is organized as follows. Section [II] describes 
the computational model of myelinated nerve fibers, the 
electrical potential field generated by the stimulation, and 
the underlying afferent activity traveling in these nerve 
fibers. Section [III| reports the results of our simulations: We 
identify different types of interactions between stimulation 
and sensory activities and we define a measure of selective 
relay of sensory neural activities. Section IV discusses the 
physiological interpretation of these results. 
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Fig. 2. Our model of extracellular electrical stimulation of myelinated 
nerve fibers in the dorsal column includes the underlying sensory activity 
as a current source at one end of the nerve fiber. 


II. Model description 

In this section, we describe our simulation test bed of ex¬ 
tracellular electrical stimulation on myelinated nerve fibers in 
dorsal column with underlying sensory activity (see Fig. [2]). 
The model is fairly simple: this is both an asset (in the math¬ 
ematical analysis) and a limitation (in the detailed modeling 
of more complex structures). As our main interest is to study 
the effect of interactions between neurostimulation-evoked 
and underlying activities, this model provides a starting point 
toward the elucidation of realistic mechanisms. 

A. Myelinated nerve fiber 

A myelinated nerve fiber is a cylindrical active membrane 
(axon), tightly wrapped in an insulating myelin sheath. This 
myelin sheath is interrupted periodically, leaving short gaps 
where the axonal membrane is exposed. 

Following McNeal’s model [6], a myelinated nerve fiber is 
represented by an (infinite) series of compartment elements 
linked by intracellular conductances. The dynamics of the 
membrane potential V- L = V- ni — V? xt at node i (where V- nt 
and V? xt are the intracellular and extracellular potentials) 
read as follows 

c m Vi + y2 kk = a (v-i - 2 v+ v i+1 ) 

keJC 

+ Ga (V™\-2Vr x +V%\) , 

where G m is the membrane capacitance and G a is the 
internodal conductance. Ionic currents at node i include 
a sodium, a fast potassium, and a slow potassium ion 
channel, as well as a leakage current across the membrane 
based on the Frankenhaeuser-Huxley model [7], adjusted to 
experimental data of human sensory fibers at 37 °C [8]. A 


complete description of the fiber model and its parameters 
is presented in [8] and [9]. 

To numerically compute the response of a finite fiber, it is 
usually assumed that no intracellular axial current flows at 
the end nodes (‘sealed-end’ boundary condition) [10]. 


B. Electrical potential field generated by the stimulation 


The extracellular medium surrounding a nerve fiber is 
composed of different regions of the spinal cord (epidural 
fat, cerebrospinal fluid, white matter, grey matter, etc.), 
which have different conduction properties [11]. In addition, 
the electrode can also take various shapes (single contact, 
array of contact, etc.) and various configurations (monopolar, 
bipolar, etc.) [12]. 

However, the extracellular medium may be assumed to 
be infinite and isotropic with the electrode represented by 
point sources at the center x c - of each contact. Therefore, the 
electrical potential field at time t and position x is given by 
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where /j n is the current of point source j and p m is the 
extracellular medium resistivity. The extracellular potential 
at node i is given by V? xt (t) = p(t,Xi), where Xi is the 
position of node i. 

The stimulation current input I st1 (t) consists of the repeti¬ 
tion, at a constant frequency, of symmetrical biphasic pulses 
with an amplitude of 2.5 mA and a duration of 350 ps. We 
consider stimulation frequencies ranging from 0 to 250 Hz. 


C. Underlying sensory activity 

The dorsal column contains nerve fibers that relay pe¬ 
ripheral sensory inputs to supraspinal centers. In normal 
conditions, these myelinated nerve fibers that originate in 
the low-threshold primary sensory neurons that mostly signal 
non-noxious sensory stimuli: proprioception from skeletal 
muscles and mechanoreception from the skin. However, 
in pathological conditions, mechanical hypersensitivity after 
injury may also be signaled by abnormal activity in dorsal 
column fibers [13], [14]. Therefore, the spiking activity in 
these fibers spans a broad frequency range and exhibits var¬ 
ious patterns [15]: regular spike discharge, regular discharge 
of doublet spikes, bursting patterns, sporadic activity with no 
regular or predictable firing pattern, etc. 

The presence of underlying sensory activity in the nerve 
fiber is represented by replacing a ‘sealed-end’ boundary 
condition by a current source at one end of the nerve fiber. 
Therefore, the dynamics of the first node becomes 

c m v 1 + y2kk = G a (v 2 -v 1 ) 

keK. 

+ G a (y 2 ext - vr l ) + J sen (t), 

where the input J sen (t) represents the underlying activity. 

As a first step, the underlying sensory activity input J sen (t) 
is modeled as a Poisson train of square pulses with an 
amplitude of 5nA and a duration of 1ms. Therefore, the 
instantaneous firing rate A sen is assumed constant, ranging 
from 0 to 100 Hz. 
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Fig. 3. Different interactions between stimulation-evoked activity and sensory activity are illustrated by their scatter plots (horizontal axis is time, vertical 
axis is position along the nerve fiber). Each dot corresponds to an action potential at a given time t and position z along the fiber. Red and green dots 
indicate the fiber response to sensory and stimulation inputs, respectively. Sensory action potential waves travel orthodromically from the bottom to the 
top of the fiber; stimulation-evoked action potential waves travel orthodromically and antidromically from the center to the extremities of the fiber. Gray 
dots indicate the fiber response in the absence of interactions (collision or loss of excitability), that is, the action potential wave that would be produced 
by the corresponding input if it was not perturbed by the activity induced by another input. 


III. Results 

In this section, we show results of our simulation test 
bed for a monopolar electrode placed at 3.5 mm from the 
center of a 10 cm-long, 10 pm-large (diameter) fiber. All the 
numerical simulations and analyses were performed with 
MATLAB, MathWorks. 

A. Interactions between SCS-evoked and sensory activities 

We identified different interaction types occurring between 
simulation-evoked activity and underlying sensory activity. 
In Fig. [3j each panel represents a typical scatter plot for 
each interaction type. In these panels, a dot is an Action 
Potential (AP) at node position z along the fiber and at 
time t. The stimulation input triggers the orthodromic and 
antidromic propagations of an AP wave (green dots) from the 
stimulation position at the fiber center toward the fiber ends. 
The sensory input triggers the orthodromic propagation of an 
AP wave (red dots) from the bottom to the top of the fiber. 
Gray dots indicate the fiber response of each input in the 
absence of the other, that is, artificially without interaction. 

The interaction type depends on the timing of both inputs 
generating these activities. 

• A collision occurs when the orthodromic sensory AP 
wave and the antidromic stimulation-evoked AP wave 
meet and cancel each other. It happens if a sensory pulse 
is triggered slightly before or after a stimulation pulse, 
that is, tf n E [tf - A t™\ tf + At™ 1 ). 

• A loss of excitability sen-sti occurs when the stimula¬ 
tion input doesn’t excite the nerve fiber due to the recent 
passage of the orthodromic sensory AP wave. It happens 
if a sensory pulse is triggered before a stimulation pulse, 
that is, tf n E [tf - At™ 1 - A t\ tf - At™ 1 ). 

• A loss of excitability sti-sen occurs when the sensory 
input doesn’t excite the nerve fiber due to the recent 
passage of the antidromic stimulation-evoked AP wave. 
It happens if a sensory pulse is triggered after a stimula¬ 
tion pulse, that is, tf n E [tf + At c _f, tf + A tf 1 + At 11 ). 


In addition to these interactions, we also identify ‘self¬ 
interactions’, that is, interactions between activities generated 
by the same input. 

• A loss of excitability sti-sti occurs when the stimulation 
input doesn’t excite the nerve fiber due to the recent 
stimulation of the fiber. It happens if the frequency of 
stimulation is too high, that is, tf +1 E [tf + At m ). 

• A loss of excitability sen-sen occurs when a sensory 
input doesn’t excite the nerve fiber due to a recent 
sensory input. It happens if two consecutive Poisson 
pulses happen too quickly, that is, tff t E [£^ en + At w ). 

B. Selective relay of sensory neural activities 

Depending on the physiological origin of the underly¬ 
ing sensory activity, we may want to modulate differently 
the relay of this input with the electrical neurostimulation. 
Typically, we want to block pathological sensory activ¬ 
ity (as mechanical hypersensitivity) but we also want to 
keep relaying normal sensory activity (as proprioception or 
mechanoreception). Let us defined reliability as follows 

. number of relayed sensory pulses 

reliability =----- — r - — -, 

total number of sensory pulses 

where relayed sensory inputs are sensory inputs that travel 
from the bottom up to the top of the fiber. 

Figure [4] shows the general decrease in reliability of 
sensory pulses as a function of the stimulation frequency 
for different instantaneous firing rates of sensory input. In 
addition, the simulations show a rebound (localized increase) 
in reliability around 150 Hz and then a decrease again. We 
hypothesis that this rebound is due to ‘loss of excitability sti- 
sti’. Indeed, 150 Hz corresponds to the frequency at which 
the nerve fiber doesn’t respond anymore to each stimulation 
pulse due to a too high stimulation frequency. In addition, 
Figure [4] shows that a higher reliability is achieved if the 
sensory input has low instantaneous firing rate. 

Remark 1: The results presented in this subsection are 
drawn from three simulations of this model with a stochastic 
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Fig. 4. The reliability of sensory input decreases as a function of 
the stimulation frequency. Dots represents the reliability for 3 different 
realizations of the Poisson input for each stimulation frequency. The line 
represents the mean of these 3 reliabilities. 

input (the sensory activity is Poisson). These results have to 
be confirmed by running Monte Carlo simulations, that is, 
running the same simulation multiple times for different re¬ 
alization of the stochastic input and computing the expected 
value of the reliability. 

IV. Discussion 

Dorsal column is the primary target of SCS [16]—[18]. 
The dorsal column contains axons that originate in the low- 
threshold primary sensory neurons that mostly signal non- 
noxious stimuli. These myelinated axonal fibers relay periph¬ 
eral sensory inputs to supraspinal centers and have collateral 
branches projecting to the dorsal horn. Since mechanical 
hypersensitivity after injury may be signaled by abnormal 
activity in A/3-fibers [13], [14] and since the dorsal column 
contains axons that originate in these neurons, inhibition of 
A/3-fiber inputs may partially contribute to SCS analgesia, 
especially for inhibition of mechanical hypersensitivity. For 
example, SCS may reduce activities in dorsal column fibers 
from reaching second order neurons in brainstem, including 
cells in the gracile nucleus and cuneate nucleus. Thus, SCS 
may interfere and alter the information coded by physi¬ 
ological sensory inputs. Or instead, AP interactions with 
afferent inputs may also change the pattern of SCS (for 
example, frequency) and inhibit antidromic APs to activate 
the segmental spinal network for pain inhibition. 

Our simulation test bed is a first step toward a better 
understanding of the effect of spinal cord stimulation in 
relay of sensory input in dorsal column fibers. In particular, 
we investigated the impact of the stimulation frequency on 
sensory input with different instantaneous firing rate (fre¬ 
quency content). We identified different types of interactions: 
collision and losses of excitability. In addition, as expected, 
a higher stimulation frequency leads to a low reliability of 
the sensory input. The choice of an optimal stimulation 
frequency may results from the dual objective to relay a 
normal sensory input and to block a pathological one, where 
the instantaneous firing rates of each input are different. 

In the future, we plan (i) to consider more complex sensory 
inputs than Poisson, such as doublets or bursts, (ii) to derive 


an analytical expression for the fiber reliability, and (iii) to 
augment the dorsal column model to include collateral fibers 
to dorsal horn and dorsal horn circuit itself. 
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